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The motion of a wingless rocket in à vacuu.v dnovt 2 spyerice! 
non-rotating earth describes the cliiptic orbit of د‎ .vwaterial point of 
masg in a central field of force. The corrections for the rotation of 
the earth are applied to this solution to determine the effects n op- 
timum range. Since there is no slinple mathernatical s slation for 
optimising the range, it is T to obtain ronpgeo for several 
angles of elevation at each initial velocity and then find ос орёл 
range by using interpolation formuias. ‘my е сасе ої {ірі at 
the equator is computed although the formulas and the com putationa: 
procedure outlined can oe applied at any latitude. 

Firing to the Last cives the greatest range for a given initics 
velocity and firing to tae West cives the least range. The angle of 

elevation is lowest firing Mast and highest firin; West, ulfvoush for 
iow velocities the difference is not gre at and for the 5,202 uisb. وق‎ 
ivtion & constant angle of clevation gives ranpes varying, at net omy 
nne-tenth of a raile froin the maximus: for any direction of fire. 


ever, for the 15, 000 rnph. solution all the factors are critical in ор 


taining the maximum range. 
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OYMBOUS AND DIA 12 
Ф = radius ol the spherica! сагу, 3980. ® OC 
M. z initial velocity relative to the cort! 
Vo = tow speed of the earth 
م8‎ = non-rotationa!l earth sravititional acceleration 
2 
= р R 
۸م‎ Bo 
دس‎ = angular velocity sf the earth 
7 = spatial velocity 
Wy = relative angle of clevation 


spatial angle of elevation 


F 
! 


relative angle of fire, :neasured соотбексіосісуіге 
from Баве 


^S 
~ 
T 


6 » = spatial angle of fire, measured clocoviise fvoc F ut. 
8 و‎ = relative bearing of point of fali, rceasured countcre- 
clockwise fron: rast 


semi-major axis of the elliptic trajectory 


p 
0 


cC 
i 


6emil-eminor astis of the elliptic trajectory 


3 = eccentricity of the elliptic trajectory 

n = one-half of the angle subtended at the center of the non- 
rotating cartn by the trajectory 

te = time of flight 

a z co-latitude of the point of fire 


longitude of point of fire, arvitrarily chosen ac ۸ 


= latitude of point of fall 


or > 


= longitude of point of fajl on non-rotatin§ cart" reintive t3 À 
(9 also used as an angle at the center of tuc curt) uf ۔ وگ‎ 
fined in Fig. 1) 


longitude of point of fall on rotating cart relative t A‏ = ما لها وا 








س = 


= 9 -Wt. - > longitude of point of fail эт то б M 
relative t5 a point of fire having Гетеа 2 

























- Ф || K ° == 
ve,“ ‘ime of F lig at, and Point of Dali 122% 


Various Angles ol Elevetion at v Lal , 090 mph. 





for ۱ 
= E 
Range, Time of Flight, and Point of Fall Deerin, 
for Various Angles of Llevoation at v s 10,9000 ipi. 


Range, Time of Flight, and Point of «all fearing 
for Various Angles of Elevation at v = 15,000 mph. 


Rotating “arth Optimum Range, Time of Flight 
and Point of fall Bearing for v=5000, 10,000, and 
15,090 raph. 


Non-rotating Marth Solution for Optimum Mange 
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Diagram on pag 
Diagram on page 6 
Diagram on page & 
Range vs. Angle of Elevation, v = 5,000 mph. 


10, 909 зарћ. 


u 


Range vs. Angle of Ulevation, v 
Range va. Angle of “levation, v s 15,000 rape. 


Time of Flight vs. Angle of “levation, v = 5,0602 soph. 


0 


Time of Flight vs. Angle of Elevation, v = 13,990 mph. 
Time of Flight ve. Angle of Elevatinn, v = 15,0060 mph. 


Deviation of Point of Fall from Firing Line for 
v = 5,000 mph. 


Deviation vs. Firing Line for v = 5,000 mph. 


D «viation of Point of Fall from Firing Line for 
у = 10,900 mph. 


Deviation vs. Firing Line for v г 10,029 тар. 


Deviation of Point of Fall from Firing Line for 
v = 15,009 mph. 


Deviation vs. firing Line for v = 15,000 mph. 


Optimum Range and Angle of Mlevation vs. Direction 
of fire for v = 5,000 mph. 


Optimum Range and Angle of Flevation vs. Direction 
of tire for v = 10,000 mph. 


’ptimum Range and Angle of Dlevation vs. Direction 
of Fire for v = 15, 000 mph. 


Optimum Range vs. Angle of Elevation aud initial 
Velocity 


zine of flight vs. Optimum Range 


Deviation vs. Initial Velocity 
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іе motion of & winglees rocket in a vacuul. dso0ut 2 spoer- 
ical non-rotating earta describes the ciliptic orbit of à v. aterisi 
point of maso ina central field of force. The equations for Buch a 
trajectory are well defined and can be reacily optimized, out wien 
considering the effects of a rotating earth there is no civ plc плав = 
ematical eolution for optimizing the range. 

To analyze these effects i¢ is necestary tu solve Mé equas 
tions for various initial velocities at various angles of elevation 
and optimize the results by interpolation formulas. 

The difficulty of carryin} put toore than eighty consecutive 
and, generally, dependent co:nputations and arriving with the re- 
quired accuracy in the answer is inherent in thé provlems. The 
computational procedure outlined in Part Y and carried out in Lable 
I ior an initial velocity of 5,000 mph, firing Northeast at the eqha- 
tor is the end result of izany attempts to shorten the co:rputations 
yet preserve the accuracy. The procedure can be carried out for 
any latitude of the point of fire and is not restricted to tie equator 


casc. 








PART ii 


MOTION OF A MAYERI cL PORNT OF ао 
Т] А. Gra МАГ, | + Tet шы C LN IRC a 


Let 0 he the center 












of the earth having a radius 
R, then the attractive force 
fnr masg m is: 


- 


Arm 
m w سب‎ 


12 


where {1 is a constant and r 





is the distance from 0 to tne 


macs. rig. i 


Г = Е = 
Р.Ж. < -( — {-dr} = < ==, if the potential eneryy at 
oF 
infinity is chosen as zero. 
The inetic ول‎ is: 
K. E. = 2 ~ where ^ ^ e (em) + (2293 | 
I 
: 2 4I ک‎ 
cm she [m سر‎ 
In a conservative field the total cnergy is constunt, hence: 
| _ 2 № 1. Га a ۳ | 
omy - ІН 5 me de + im. (1) 
^ 


A further condition ib that tne moment of r5o'nentars reécain | 


constant about 0, hence: 


ха (r E) r = r. (v. coo V) R (2) | 
| 


سا < = = 


* Reference (1), page 35 








n 51 л ———" 3 з жаг а 
e (1) and 12) Geécribe the motiaa ai sl. 


I _ ar d$ ж ы 
uttins Toe 2a 22 in Ka. (2) 


T^ 
мн cas V 





de = dr I | N. | 
dt 2 de , and using this relation 


im Eq. (1). 


» 
Нето бое лу En + R^ vf cos" Y ele _ 8 о 5ه‎ 











€ r 2 e. 
> E 
EU. 2 1 
7, === 
ау „2 ې‎ E pe 
т Av. cos V - 
- г % G ۲ e 
3 E Avv COD V 
> 
Rv cos WV 1 





"те су C‏ - رلوب 
P У =‏ 





سه رع 


2. Rv 5608 ү | 2 
^ "OX ۴ СА 7 яё ! ا‎ ҢҮ گآ حا‎ 
WA J 





Integration gives: 





TK 
Ry cosW : : 
7 Ху- cos Y 


cos (8 + 8°) = : 2; 1 
۳ - ٦ اسم‎ 2 


R ۷ 625 Y 


t 


wnere Ө! іс the integration constant. 
Solving {or r, 


Мм, СОВ y 
>,» 2 21" E t 
Ee Goan de ~ e2 Ag ++ 7 Cot (oy is ) 
Vv cosy Q ЈА + - 2 
2 +. Y Сов Y 
a 
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^1! 


For r to have its maximum value Q = 0 and therefore 6 


must be y giving 


Rv cos Y 
ی‎ M "3 к” - 
"v cos vU < admo ۳ A ہد‎ = 
i Mo SN о К Rev cub ү 


Now the general forrn for the equation of a conic section i8: 


rs РЕ » where p and € are constants. 


Let 8, be the gravitational constant for a non-rotating earth 


at the surface, and putting K = د‎ then 


Bm 1 
ua و رس سی‎ 2 » therefore 
2 
Уо сов ү 
бо 








| Y. 2 ë 2g, R 2 
l~ ‘Ree cos W (1 - "+ + 1 сов 0 
О 


ihe eccentricity is then 


2g X 2‏ و 
;7% 7 
[e eye S |‏ 
У‏ 9 
J‏ 


and if € < 1, tbe orbit is an ellipse 
if E£ zl, the orbit is a parabola 


and if £»1, tbe orbit ie a hyperbola. 





lor € ts be less than 1, then 





eg," 
Ë : = 
= Гот у < ری‎ 
У 
О 
7 = ڑا‎ then € = sin WW, ana if firing horizontally 


(М = 0) the eccentricity vanishes giving a circular orbit about the 
earth. 

Ail uf the following computations are for velocities less 
than the above "escape velocity", thus only the elliptic orbit is 
considered. 

Using the polar forrn of the elliptic equation®, the half 


major axis, a, is: 


ws‏ سر 
м‏ “~ 
E E 2‏ 82 = 2 
The half minor axis, b, is:‏ 
“ 1 34 3 


у.“ 2g R > 
د و‎ 


"15 э | 
ay R COS WA > 
JA 





or 8 
„2 К NÖ we 2 
AA 


= => че се <> 


#* Reference (2), Chap. Vi, §17,17. 


GLLIPTIO TRAJECTORY ABOUT 








Е _ 
= 


Jin & vacuums about a spheri 
Г —— 


та Ww 
cal, non-rotating carth. 
— о. 





The foci are O anî ان‎ 
where O igs the center 


of the eartr. 
т 
Y = > = W 
а =т= 2у = 2 V 


O's 2h - AS 


n! | = 2га - Ж = | 


/ 
| 
| | 
I 454 
| 
| | 
| 
| | 


Fig. 2 


5 1 
“ence: 0 = ES + (4 „ 21, ] cos 2 y] 2 


: Dein Ё 

[87 + 17-28 Дсова У 

Phe arc lengt:. AB is the range вес ۶ 
Range = сир 


The time of flight to traverse the distance from A to B is 


given by*: 


"P a ФӘ ар Фә 


"i.eference (23, Chap. VI, 412. 
i 




















ART IV 
TE20T YU ROTA TIM DARTH 
“з determine the effect of the rotating carth, the latitude 


and longitude of the point of fall on the non-rotating earth ic first 


ound to bet (вее Fig. 3): 


Z 


7 


NA 


X 


Latitude; 


DEN! ۱ = 
ф = sin “(cory cosa - siny cos Û i: ina) 


i.ongitude: 
Ба s siny ein ß 


oin2ysin®ß +Hcinycos 6 сова -cosysino)" 
The preat circular arc AP is the projection of the trajec- 
tory on the non-rotating earth. Using the equations derived above 
with the initial velocity as the vector sum of the relative initial vce- 
locity with respect to the earth's surface and the tow sneed due to 


the rotation of the earth yields a trajectory wits respect to & non- 


rotating syste: of coordinatef with the earth's center 26 the origin. 


* Reference (2), Chapt. VI, 920 








۳9 





Since t, time units are required to go from A to E then if ته‎ ۸+08 


ke angular velocity of the earth's rotation, the ear thi rotates about 














А 
DANT V 
COMPUTATIONAL PROCEDURE FOR ROTATING ZANT 
Since the radius of the earth R appears in many of the egua- 
tions and the velocities are large, the calculations are considerably 
simplified by taking velocity units as radii of the earth per hour. 
Thus any place that R appears in the above equations unity replaces 
it. У | 
The initial values and constants used are: 
Yor 5,000 mph. relative initial velocity, 
v, = 1.26298027 Radii/hour 
For 10,000 mph. 
ма = 2.52596054 Radii/hour 
For 15, 000 mph. 
Ve = 3. 78594082 Radii/hour 
The tow speec of the earth at the ecuator is: 
· У, 59. 26251645  Radii/hour 
The non-rotational gravitational acceleration at zero height 


is taken as that at 45° latitude on the actual earth and then corrected 


for non-rotation. 
и, = 19.99727) Radii/hour” 


and 
2 2 59 2. 
М =g F =19.997270 Кади /hour 
А = 4. 4718307 дайн? /2 /hour 


The angular velocity of the earth is taken as 


© = .26251645 Radians/hour 


Colusans 2 through 22 of Table 1 give the vector surmmation of 








89 
relative initial velocity and tow speed of the earth. The spatial ves 
amana... 
the spatial angle of elevation, W ىد‎ si 4 the spatial angle 


locity, v 


o! 
of fire,  ,, arc used in the remaining calculations. (Note that f, 


is measured counter-clockwise from East, and 6, is measured 
clockwise from South solely because oi computationai ease.) The 


equation for a using X = l becomes: 


1 
а = 2.086 
A 
and is computed in Colurangs 23 through 25. Columns 26 and 27 are 


MED 


required later on and give а2 and a 
The equation for J solved for in Column 28 is: 
Д = 228 + 1 
š ) 2 
Column 29 gives f 


The equation for نا‎ solved for in Column 30 is: 
5^ + =. 2 (2 ) 
۱ = асов V. № 
The equation for £ “° 2 solved for in Column 31 and 32 is: 

21 pl 
£ = исту 

a 

Related functions of € are given in Colurnns 33 through 36. 


The equation for X) solved for in Columns 37 through 43 is: 


ia we T cin 2 Vg 


2 
1 + 1-2 f coe алу. 


The ec p for t, solved for in Colum ns 44 through 54 is: 


2 


/2 
= ^^ | 6. 2831853 - A tan” 1-5. tan (1. 5707062 - 5) 


-1 


e 


+ 2E sin | 2 tan (1. 5707963 - | 











lien 
Columns 55 through 85 are spherical trigonometric formulas 
and are self-cxplanatory. (itote that a is the co-latitude of the point 
of fire). 
T.ble Il shows tne interpolation procedure for optimizing 


the range. The Gregory-Newton interpolation formula ів: 





12 ) 1 2 

yey +nay + HEL ay ې‎ E دو‎ 
and 

d 1 (2n-1) ,2 

az Қ (һу, 2529, Жа, ІР ДА е о e] 

Where Ay, is the first difference, Ay. the second differ- 
š х-х 
ence, etc; Ax is the constant difference between x's; n = > р 


26 


6 ۱ 
Hence for a = 0 {maximum range), 


Ay 
n = .5 < کہ‎ 
^ yg 
and 
охбу. 
x =X +.54X- u. 
ға 0 A^yg 


The maximum ranges and the corresponding angles of ele- 
vation thus found are shown graphically in Figs. 4, 5, 6. 

The time of flight, te; and the azimuthal angle of bearing of 
the point of fall are found graphically from Figs. 7, 3, 9, and Figs. 
10, 11, 12, 10a, ila, 12a. 

The pertinent results from all the computations are listed 
in Tablesll, IV, and V, while the maximums of these Tables are 


listed in Table VI. 
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PART Vi 





case at cach velocity it is only necessary to maximize the equation 





previously given for N : 





P ein 20۷ ےس‎ 


- Е тр aT 
ы. 1+ ۱*۰ 21 6062 با‎ 


for N max then: 


cos 2 V = ] 


and 
sin n = q 
where, 
و ود‎ не 
> d 
m beu 
^^“ 


Table VI shows the calculations and results for the three 


velocities used in the non-rotational earth solution. 
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PART VI 
DISCUSSION OF RESULTS 

Jigs. 13, 14, and 15 show cross plots of Digs. 4, 5, and 6 
as compared with the non-rotational values oí Y and range. 

ius. 16 and 17 give a comprehensive graphical solution to 
a particular firing problem, Given the range and direction of fire, 
then the angle of elevation, V , the initial velocity, VP and time 
of flight, t,, can be found. 

Fig. 18 shows the change in deviation of the point of fall 
with increasing initial velocity for the case of firing Northward at 
the equator. The left deviation at low velocities is due to the rocket 
having à constent angular momentum about the axis of the earth 
during its flight and the trajectory remains mostly outside oí a hy- 
pothetical cylinder having the earth's axis and being tangent to the 
earth at the equator. Yor the rocket to land on the same bearing 
as it was fired, the angular momentum would have to vary directly 
ав (һе rocket's distance и the earth's axis throughout ita flight. 
tence portions of the trajectory outside the cylinder cause left 
deviation while portions inside the cylinder cause right deviation. 
Due to the flatness of the earth's surface for short ranges there is 
an almost linear increase in left deviation for initial velocities up 
to 6,000 mph. As the velocity increases, the angle of elevation 
decreases so that less and less of the trajectory is outside of the 
cylinder. At an initial velocity of 11,600 mph the trajectory is such 
that the left deviation portion cancels the right deviation portion. 

At higner velocities the deviation increases rapidly to the right. 
| А similar analysis can be applied to the range effect and tu 


cases of firing at latitudes other than the equator. 
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TAMLE i* 


notetiny Darth Solution for Range at VY „=. 650, .725, ск) ГЕС, 
u = 5000 mph, Firing Northeast ( A, = v/4) 


sin(3) сов (3) 
T 


sin سس‎ 






cos 







5in(2) 
бір Pr | 


.T07106676 


1. 00543813 






. 70710970 





94534009 












.2799%ҰЫ65 .7071067С 












1. 2054454 






‚97345849 . 9476409 





. 19710967) 








. 70710679 .93097269 . 9667101 1. 1460996 


















. 88471851 1.081601? 





. 70710670 . 7821269 





. 0375073 







رر . 72 .91 ° 





, 76659204 





. 9907107 





< > = سا 


т ` 


*Peferences [5) and 14) vere oget in these computations. 








Edi. 











; 21 ! 22 
асобе 
| 42 с | Gi с; с: 
| cos Ws | sin2 Ws | cos Z Vs 
| ‚ 7132424 | ‚ 9044947 | ‚ 4264848 ‚ 1018793 1. 3981207 
' | | 
| -6518844 | .9527458 | 3037688 . 1007635 1. 8992365 | 
| | | 
| .5876608 | .9845112 !  .1753216 ‚ 0995491 1.9004509 








| 25 26 27 29 

l 

x 1/ (2 2094 | rx ШЕГЕН ef 
5 f №: | 


‚ 5268369 227750 7) | 0 ‚ 0536738 ‚ 0028809 











. 0530548 . 0028148 






£ 
. 5265274 . 2112311 | . 08543702 





| 
| 
| 


. 0523518 









. 5261909 . 2768769 . 08535513 . 0027439 














. 0382324 213792062 | . 8620738 . 9284793 . 0715207 
. 0345856 . 1247537 | . 8752463 . 9355460 . 0644540 
. 0307827 .1111783 | . 8888217 ‚ 9421734 . 0572266 





E e os [os آاںے‎ 
__-з |Ја | 8 


1. 9284793 . 19257883 . 0485477 1. 0028809 . 0457821 





1. 9355460 . 1824834 . 0505477 1.0028148 0322325 





| 
| 1. 9427734 . 1716278 , 0515705 1. 0027439 . 0183673 















43 
ыш 42 
| = 7 =. 
. 9570988 . 9783143 : 0496238 ۳ ‚0496442 | š 0248221 
. 9705820 . 0513305 . 0256653 





. 9851812 | . 0513080 


| 
| | 
| | 
.9921575 | .0519781 | .0520015 |  .0260008 ` 


A nr AA 





. 9843766 


| 45 ۱ 6 1 ат | ав 09. 
1. 5707963 ı 45 жу tan” (4) | 2 : 
En | سا‎ зы ИЕ 
7 2 | tan T/2 | | | 


| 
| 
| 
| 
| 
к 1 SE II 
1.5459742 | 40.27839 | 7.7567640 | 1. 442583 38 | 2 8851676. 
| | 
{ 











' 
— | ممت ممیت ممل مہ نے ن ی سے ف 5 Е‏ 
































| | 

1. 5451310 | 38. 95455 | 7.1085587 | 1.4310379 2. 8620758 
| 1. 5447955 | 38.45164 | 6.5993704 | 1.4204108 |  2.8408126 | 
ED (| St 54 | 
= 
| e, Om А 
u, 7 ui sin u, | 2 Esin u, | | lum t; | 
. 5128501 E maro T .9838196 | 0841287 | 
‚ 5590337 | . 2158912 | .5162175 | 1.0752515 | 0918663 
| .6015421 | .2962567 , .5586059 | 1.1601480 0990246 | 
l — 1 
sin(56 | cos (56) | 
Sin y и | Ë _ сову 7| 
.0991254 |  .9950749 | 
‚0241164 | .1026610 | (1024808 | .9947349 | 
)" | .1040030 | m/2 10381566 9945966 ` 





۰ 
Rte ate Ds A DITA و یچ‎ ee D پس -- الا‎ + o. ۰ .. ` e ...- + 


*No column 57 
















& у 


203 ф 


Се а É FT GE کے‎ DASE “47% w وٹ‎ > 
-, , 2620 و‎ 52. 4526 = ۰ . G 7 23 = د راګن ار‎ 
geo fs # ç A EN IZ ` e a р гу ہی‎ 7 K qui еа = 
د چا‎ 9% 0 ( e 14 о Ke > Y (12 e Cag 424%“ 01 e < 01 21720 


=. 2597210 , 0597210 . 0597566 . 9982.13 


-, 0333539 


- 0550599 


. ۵6 









. 0354078! . 6339768 7111697 

















. 1883852 


80595903 


. 7911466 





рар Ue ap ۳ | ` ۳ € - | » ۱ 
+ Colo, 52 through 75 must be weed ih place of Cale. nd ©! 
ы А T x : m у Р А = + ! 7 - ~ N > f 
& Г) tado 12 P accurate iy OM 7/4 <. 8. < Зт /4 ñ 
^ x D > 4 " i 















VAL Pal. 1 





Optimum Капыс Solution for Nesulite of Гаоіе I 


۷ 


z 5000 mph. В. = 7/4, 
¥ | 


=. 











> وہ‎ 
x 


TABLE IH 
Range, ¿une of light, and Yoint of Tall Bearing for various 
Angles of Klevatian at v = 59000 ۰ 


v, = 1. 262980 (5000 mph) 


4 : ۳ 
> 























la 0 . 65202 .0848425 | 328.715 
b y · #2590 ‚ 0925949 | 335.402 
C G . 90000 . 0997507 | 334. 325 











. 0990246 | 332. 847 

‚ 0824536 | 321.443 1 

. 0902379 | 329.433 1 

‚ 0973135 | 329. 324 1 
65060 . 0808777 | 316.532 2 
72509 ‚ 0355023 | 524. 321 2 
80000 . 0956602 | 325. 859 2 
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TABLE 1۷ 
Range, Time of light end Point of Yall Bearing for arius 


Angles of !')evation at v = 10,967 mpi. 


de 


7 = 2. 52596054 (10,000 rnpn. } 


r 
la 75 


٤ 
a 0 . 57500 . 2063123 


- 

















hs 





. 2209331 
. 2333689 
. 2102985 
, 2277055 
2437099 
2005300 


, 72500 . 2130675 








. 90500 . 2341679 





1968131 . 0ء ۔ 
72500 . 


. 59099 
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